Differential Elliptic Flow in 2 - 6 AGeV Au + Au Collisions: 
A New Constraint for the Nuclear Equation of State 
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Proton elliptic flow is studied as a function of impact-parameter 6, for two transverse momentum 
cuts in 2 - 6 AGeV Au + Au collisions. The elliptic flow shows an essentially linear dependence on 
b (for 1.5 < b < 8 fm) with a negative slope at 2 AGeV, a positive slope at 6 AGeV and a near zero 
slope at 4 AGeV. These dependencies serve as an important constraint for discriminating between 
various equations of state (EOS) for high density nuclear matter, and they provide important insights 
on the interplay between collision geometry and the expansion dynamics. Extensive comparisons of 
the measured and calculated differential flows provide further evidence for a softening of the EOS 
between 2 and 6 GeV/nucleon. 
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For several years now, the study of nuclear matter at 
high energy density has held the promise of providing 
valuable insights on the nuclear equation of state (EOS) 
and on the predicted phase transition to a quark-gluon 
plasma (QGP) |]-| ■ At AGS energies of - 1 - 14 AGeV, 
elliptic flow has emerged as an invaluable probe of high 
density nuclear matter This flow has been at- 

tributed to a delicate balance between (i) the ability of 
compressional pressure to effect a rapid transverse ex- 
pansion of nuclear matter and (ii) the passage time for 
removal of the shadowing of participant hadrons by the 
projectile and target spectators (|] . If the passage time is 
long compared to the expansion time, spectator nucleons 
serve to block the path of participant hadrons emitted to- 
ward the reaction plane, and nuclear matter is squeezed- 
out perpendicular to this plane giving rise to negative 



elliptic flow. For shorter passage times, the blocking of 
participant matter is significantly reduced and preferen- 
tial in-plane emission or positive elliptic flow is favored 
because the geometry of the participant region exposes 
a larger surface area in the direction of the reaction plane. 
Thus, elliptic flow is predicted and found to be negative 
for beam energies < 4 AGeV and positive for higher beam 
energies 

Recent theoretical studies of elliptic flow have sug- 
gested a sensitivity to the pressure at maximum compres- 
sion [Q|||llJ and thus to the stiffness of the EOS, and to 
possible QGP formation ||. Despite this sensitivity, the 
commonly calculated patterns for elliptic flow very often 
do not constrain the EOS uniquely. This being the case, 
it is important to investigate additional experimental ob- 
servables which may provide more stringent constraints 



1 



for the EOS. Here, we investigate the utility of differ- 
ential elliptic flow measurements Viip) and ViipjPT)-, as 
possible constraints. 

The measurements were performed at the Alternat- 
ing Gradient Synchrotron (AGS) at the Brookhaven Na- 
tional Laboratory. Beams of 197 Au at EBeam — 2, 4, 
and 6 AGeV pj were used to bombard a 197 Au target 
of thickness calculated for a 3% interaction probability. 
Typical beam intensities resulted in ~ 10 spills/min with 
~ 10 particles per spill. Charged reaction products were 
detected in the Time Projection Chamber (TPC) Q 
of the E895 experimental setup. The TPC located in 
the MPS magnet (typically at 1.0 Tesla) provided good 
acceptance and charge resolution for charged particles 
— 1 < Z < 6 at all three beam energies. However, 
a unique mass resolution for Z = 1 particles was not 
achieved for all rigidities. Data were taken in two ex- 
perimental runs with a trigger which allowed for a wide 
range of impact-parameter selections as presented below. 

Our flow analysis follows the now standard proce- 
dure juj of using the second Fourier coefficient, v 2 — 
(cos 2(f>) , to measure the elliptic flow or asymmetry of the 
proton azimuthal distributions at mid rapidity (|y C m| < 
0.1). This distribution can be expanded as 

dN 

— oc [1 + 2 vi cosU) + 2 v 2 cos(20)] , (1) 
a<p 

where <f> represents the azimuthal angle of an emitted pro- 
ton relative to the reaction plane. Near mid rapidity in 
a symmetric system V\ « 0. The reference azimuthal an- 
gle $ p iane of the reaction plane is determined using |lj| 
the vector Q, = Y^jtiMv^PTj/PTj- Here, p Tj and 
Uj represent, respectively, the transverse momentum and 
the rapidity of baryon j (Z< 2) in an event. The weight 
w{yj) is assigned the value < p x > / < pr >, where p x 
is the transverse momentum in the reaction plane 
The average < p x > is obtained from an earlier pass 
of an iterative procedure employed for each energy and 
impact-parameter selection. 

The orientation of the impact-parameter vector follows 
azimuthal symmetry about the beam axis. Therefore, the 
azimuthal distribution of the determined reaction plane 
should be uniform or flat. We have established that de- 
viations from this uniformity can be attributed to defi- 
ciencies in the acceptance of the TPC and have applied 
rapidity and multiplicity dependent corrections following 
Ref . . The corrections were applied for each of several 
impact-parameter selections at each beam energy; they 
ensure the absence of spurious elliptic flow signals which 
might result from distortions in the reaction plane distri- 
bution. The dispersion of the reaction plane < |<^i2 > /2 
was estimated for each impact-parameter b via the sub- 
event method |^5| . Suffice to say, a reasonable resolution 
was observed over the entire range of energy and b stud- 
ied. These estimates for the reaction plane dispersion 
serve as a basis for evaluating the dispersion corrections 



summarized in Table 1; these corrections have been ap- 
plied to the extracted flow values discussed below. 

The event multiplicity of identified charged particles 
Mfut was used for centrality selection. That is, several 
multiplicity bins were selected in the range from 0.4 to 1.0 
Mmax where M max is the point in the charged particle 
multiplicity distribution where the height of the distri- 
bution has fallen to half its plateau value |fl6|| . Impact- 
parameter estimates have also been made for these cen- 
trality selections, at each beam energy, via their respec- 
tive fraction of the minimum bias cross section. 

Figure |Tj shows representative distributions in the az- 
imuthal angle <f> obtained at the energies of 2, 4 and 6 
AGeV for mid-rapidity (|ycm| < 0.1) protons. The pan- 
els from left to right represent the three beam energies, 
respectively, and from top to bottom the three impact- 
parameter ranges of < 6 < 3, 4 < 6 < 6 and 7 < b < 8 
fm. For visual clarity, a px cut has been applied to the 
distributions shown for both the 4 and 6 AGeV data, as 
indicated. Within each &-range in Fig. |l|, the previously 
reported transition from negative to positive elliptic flow 
at rs 4 AGeV (?J is clearly seen. That is, the elliptic flow 
is negative at 2 AGeV, positive at 6 AGeV and essentially 
zero at 4 AGeV. An apparent increase of the anisotropy of 
the distributions with increasing b can also be discerned 
for the 2 and 6 AGeV data shown in Fig. [p. We attribute 
this trend to an interplay of the changing geometry with 
the expansion of excited participant matter as discussed 
below. 

Figure || shows the v 2 coefficients for the full px range, 
as a function of b for data (filled stars) obtained at 2, 4, 
and 6 AGeV in the three panels, respectively. These coef- 
ficients have been obtained by evaluating the (cos 20} for 
each azimuthal distribution obtained for a given impact- 
parameter at each beam energy. A correction has been 
applied to some of these coefficients to account for biases 
resulting from (i) low pr acceptance losses in the TPC 
for the 2, 4, and 6 AGeV beams, (ii) high px acceptance 
losses in the TPC for the 2 AGeV beam, and (iii) 7r + con- 
tamination of the proton sample at 4 and 6 AGeV . 
A procedure for effecting these corrections has been de- 
tailed in Ref. 0. That is, we first plotted the observed 
Fourier coefficient (cos 2(f)) vs. pr with pr thresholds 
which allowed clean particle separation (pr ~ 1 GeV/c). 
We then extracted the coefficients for the quadratic de- 
pendence of (cos2(/)') on pr- These quadratic fits are re- 
stricted by the requirement that (cos 2(f) ) = for px = 0. 
Next, we corrected the proton pt distributions for possi- 
ble high and low pt losses. A weighted average (relative 
number of protons in a p^ bin times the (cos 2(f) ) for 
that bin) was then performed to obtain (cos 2(f)) for each 
beam energy. The corrections which result from this pro- 
cedure are ~ 5% for the 4 and 6 AGeV beams and ~ 15% 
for the 2 AGeV beam. Subsequent to these evaluations, 
the V2 values were corrected for reaction plane dispersion 
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using the procedures detailed in Refs. 0,|TJ,|T^,|T8| . 

The i>2 values represented by filled stars in Fig. |^ indi- 
cate an essentially linear dependence on impact param- 
eter. The slope of this dependence is clearly negative 
and positive for the 2 and 6 AGeV data, respectively. 
By contrast, an essentially flat dependence is observed 
for the 4 AGeV data suggesting that the beam energy 
at which the elliptic flow changes sign is not very sen- 
sitive to b for < b < 8fm. The approximately linear 
dependence exhibited by the data can be understood in 
terms of the collision geometry and the development of 
transverse expansion within the participant matter. At 2 
AGeV an expansion perpendicular to the reaction plane 
develops over the characteristic time of d/c s while the 
spectators are present. Here, c s = y dp/de represents 
the speed of sound for a given pressure p, and energy 
density e, and d is the perpendicular distance from the 
center of the participant region to the surface. The spec- 
tator passage time (estimated in sharp cut-off geometry) 
first increases and then remains essentially constant as b 
increases over the range of interest. On the other hand, 
the expansion time decreases with increasing b due to a 
decrease in d. It is this decrease in the expansion time 
coupled with an essentially constant passage time, which 
provides the driving force for more matter to escape the 
interaction region as b is increased i.e. an increase in 
"squeeze-out" with b. The magnitude of the "squeeze- 
out" follows an approximately linear dependence because 
d is roughly proportional to 1 jb for the Au + Au impact 
parameter range 1-8 fm. 

At 6 AGeV the spectator passage time is very short 
compared to the expansion time and preferential in- 
plane emission dominates. In this case, however, the 
linear increase of vi with increasing impact parameter 
is strongly influenced by the initial spatial asymmetry 
of the nuclei overlap region or participant matter. This 
asymmetry is commonly characterized in terms of the 
width L x and height L y of the overlapping region via 
a s = (L y — L x )/(L y + L x ) ^ and can be shown to 
be nearly linearly proportional to the impact parame- 
ter for medium b values. The essentially flat dependence 
of V2 observed at 4 AGeV suggests that, at the tran- 
sition energy, the reduction in the expansion time [in 
competition with the spectator passage time] with in- 
creasing b, is compensated for by the (later) increased 
in-plane-emission from the preserved initial spatial asym- 
metry. 

The solid circles, open squares and solid triangles 
shown in Figs. 0, represent results from calculations 
with a recent version of the Boltzmann Equation Model 
BEM || which assumes a soft (K = 210 MeV), a stiff 
(A' = 380 MeV) and an intermediate (K = 300 MeV) 
EOS respectively. The calculations include momentum 
dependent forces |l{| . A comparison of the calculated V2 
values indicate sizeable differences between the predic- 
tions for a stiff and a soft EOS for all three beam ener- 



gies. For both 2 and 4 AGeV this distinction increases 
with increasing impact parameter indicating that the im- 
pact parameter dependence of elliptic flow lends a new 
and important constraint for the EOS. At 2 AGeV, the 
v-i values for the stiff EOS show good agreement, both 
in magnitude and trend, with the experimental data. At 
4 AGeV the measured values lie between the calcu- 
lated result for a stiff and a soft EOS, and appear to be 
in better overall agreement with an intermediate form of 
the EOS. At 6 AGeV the data is less compatible with a 
stiff EOS, but does not allow a clear distinction between 
the soft and the intermediate (K = 300 MeV) EOS. The 
results of these comparisons can be taken as being sug- 
gestive of a softening of the EOS as previously reported in 
Rcf. 0. However, it is interesting to investigate whether 
or not the differential flow measurement vi^b^pT) pro- 
vides further constraints. This line is pursued below. 



Figs. H compares experimental (stars) and calculated 
(circles, triangles and squares) differential elliptic flow 
V2(b,pT) for 2, 4 and 6 AGeV as indicated. At each beam 
energy, the BEM calculations have been carried out for 
the same pr and b selections applied to the data. Fig. || 
indicates good agreement between the data and the cal- 
culated results for a stiff EOS at 2 AGeV. At 4 AGeV 
the data again shows better overall agreement with the 
intermediate and soft EOS. At 6 AGeV the comparison 
indicates quite good agreement (both in magnitude and 
trend) between the data and the results from the calcu- 
lations which assume a soft EOS. The latter agreement 
is in contrast to the results obtained from the compari- 
son made in Fig. |2| and clearly indicates that the differ- 
ential flow V2(b,px), does indeed provide an additional 
constraint for making a relatively clear distinction be- 
tween the different EOS's at 6 AGeV. 



To summarize, we have studied differential proton el- 
liptic flow in 2 - 6 AGeV Au + Au collisions. The elliptic 
flow shows an essentially linear dependence on b, in the 
range 1.5 < b < 8 fm, with a negative slope at 2 AGeV, 
approximately zero slope at 4 AGeV, and a positive slope 
at 6 AGeV. This dependence can be understood in terms 
of (a) the relationship between the collision geometry, 
(b) the relative magnitude of the time for development 
of the transverse expansion, and (c) the passage time for 
removal of the shadowing of participant hadrons by the 
projectile and target spectators. Detailed comparisons 
between the measured differential elliptic flow V2{b,pr), 
and V2(b), and the results obtained from a relativistic 
Boltzmann-equation calculation not only suggest a soft- 
ening of the EOS but also indicate that differential flow 
measurements provide very important constraints for the 
determination of the EOS of high density nuclear matter. 
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FIG. 1. Measured azimuthal distributions for Au + Au 
collisions. Distributions are shown for the impact-parameter 
ranges of < b < 3 fm, 4 < b < 6 fm and 7 < b < 8fm and 
the beam energies of 2 (a, b, c), 4 (d, e, f) and 6 (g, h, i) 
AGeV, as indicated. The solid lines are drawn to guide the 
eye. 





Dispersion Correction Factor 


b range (fm) 


2 AGeV 


4 AGeV 


6 AGeV 


< b < 3 


1.71 


2.64 


4.65 


4 < b < 6 


1.22 


1.59 


2.47 


7 < b < 8 


1.26 


1.99 


2.86 



Table 1: Correction factors for reaction plane disper- 
sion for several impact parameter ranges for the 2, 4 and 
6 AGeV beam energies. 
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FIG. 2. V2 as a function of b (p T > 0) for 2 (a), 4 (b) 
and 6 (c) AGeV Au + Au collisions. Experimental values are 
indicated by the filled stars. The open squares, full circles 
and solid triangles represent V2 values from BEM calculations 
with a stiff (K = 380 MeV), a soft [K = 210 MeV) 
and an intermediate (K — 300 MeV) momentum-dependent 
EOS respectively. The identified charged particle multiplic- 
ity Mfm, is also indicated for several values of b. The solid, 
dotted and dashed-dotted lines serve to guide the eye only. 



FIG. 3. Same as Figsj2| except that a pr > 0.7GeV cut has 
been applied on the data and calculations at 2 and 4 AGeV 
and a pt > l.OGeV cut has been applied at 6 AGeV. 
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